Dynamical instability in kicked Bose-Einstein condensates 
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Bose-Einstein condensates subject to short pulses ('kicks') from standing waves of light represent 
a nonlinear analogue of the well-known chaos paradigm, the quantum kicked rotor. Previous stud- 
ies of the onset of dynamical instability (ie exponential proliferation of non-condensate particles) 
suggested that the transition to instability might be associated with a transition to chaos. Here we 
conclude instead that instability is due to resonant driving of Bogoliubov modes. We investigate the 
Bogoliubov spectrum for both the quantum kicked rotor (QKR) and a variant, the double kicked 
rotor (QKR-2). We present an analytical model, valid in the limit of weak impulses which correctly 
gives the scaling properties of the resonances and yields good agreement with mean-field numerics. 

PACS numbers: 03.75.Lm,05.45.-a,03.65.Ta,03.75.-b 



INTRODUCTION 

The production of Bose-Einstein condensates (BECs) 
in dilute atomic gases has opened up a new domain for re- 
search in quantum dynamics, since BECs are intrinsically 
phase-coherent and can be controlled experimentally to 
an extremely high degree of precision [lj . An increasingly 
interesting aspect of the dynamics of BECs is that they 
represent a new arena for investigation of the interaction 
between nonlinearity and quantum dynamics, including 
quantum chaos [2|, y, 14, la, Ifl III) B l2( • 

A BEC subject to periodic short pulses, or kicks, from 
standing waves of light represents a nonlinear general- 
ization of the well-known chaos paradigm, the quantum 
kicked rotor (QKR). The QKR has been realized using 
(non-condensed) cold atoms, permitting experimental in- 
vestigation of a range of interesting chaos phenomena 
[lfj| . The regime where the kick-period T is a rational 
multiple of 7r has also proved of particular interest: sev- 
eral studies have inv estig ated the dynamics here with or 
without linearity [1, IllL 12, [HI]. A number of experi- 
mental studies have also investigated kicked BECs 14j . 
Ensuring dynamical stability of the condensate is thus 
very important in studies of its coherent dynamics: if 
the condensate is dynamically unstable, numbers of non- 
condensate particles grow exponentially. If it is stable, 
they grow more slowly (polynomially). More broadly, the 
study of diffferent types of instability in static [15| and 
driven BECs [16| is of much current interest. 

Previous work on kicked systems^, 0, @] considered 
the onset of dynamical instability and investigated the 
relation with classical chaos. In [3(, the possibility that 
instability is related to chaos in the one-body limit was 
investigated for the Kicked Harmonic Oscillator. In 
the correlation between chaos in the mean-field dynam- 
ics, rather, and the onset of dynamical instability, was 
investigated. An "instability border" , determined by the 
kick strength K and the nonlinearity g was mapped out; 



it was then found 0] that the parameter ranges for this 
border corresponds closely to a transition from regular 
to chaotic motion, of an effective classical Hamiltonian 
derived from the mean-field dynamics. Hence, present 
understanding of onset of dynamical instability in kicked 
BECs suggests that it may somehow be related to a tran- 
sition to chaos. 

In this work, we conclude that a quite different mecha- 
nism is primarily responsible for dynamical instability in 
the QKR-BEC. Our key finding is that it is the strong res- 
onant driving of certain condensate modes by the kicking, 
which triggers loss of stability of the condensate. This 
mechanism is unrelated to the transition to chaos, but is 
rather an example of parametric resonance. In another 
context, the relationship between parametric resonance 
and dynamical instability of a BEC in a time-modulated 
trap is a topic of much current theoretical li|, 18| and ex- 
perimental interest [l9[. But to date, "Bogoliubov spec- 
troscopy" in the analogous time-dependent system, the 5- 
kicked BEC, has not been investigated. Our study shows 
that the temporally kicked BECs open up many new pos- 
sibilities in this arena. 

We find that in general, for the kicked-BEC, there is 
no single stability border: typically, for moderate K, the 
condensate restabilizes just above the stability border. 
For small K and g the number of non-condensed atoms 
N ex {t) grows exponentially only very close to a few, iso- 
lated resonance peaks. With increasing K and g, the 
number of resonances which can be strongly excited by 
the kicking proliferates and overlaps. Our calculations 
show this is associated with generalized exponential in- 
stability; however this regime is, to a large degree, beyond 
the scope of our methods. For lower K and g, though, 
we introduce a simple perturbative model which provides 
the approximate position and width of the important res- 
onances for both rational and irrational T. 

A key finding is that, for the integer values of T / tt = m 
(where m is integer) values, the focus of the study in 
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6], the onset of instability can occur at nonlinearities 
much lower than those required to resonantly excite even 
the very lowest collective mode - a key reason why the 
mechanism of parametric resonance may so far been over- 
looked in respect of destabilization of kicked BECs. Our 
model demonstrates that for this case, resonant excita- 
tion involves two excited modes in addition to the initial 
ground state mode. Hence we can explain the position 
of the critical stability border found in @, [13] ■ 

We investigate both the usual QKR-BEC as well as 
a simple modification, obtained by applying a series of 
pairs of closely-spaced opposing kicks (the QKR2-BEC). 
This modifies substantially the relative strengths of the 
resonances, and provides the added novelty that the low- 
est modes are excited by an effective imaginary kick- 
strength. It is closely related to the double-kicked quan- 
tum rotor, investigated in cold atoms experiments and 
theory [20J. We introduce a simple analytical model 
based on the properties of the unperturbed condensate, 
which gives the distinctive properties and scaling behav- 
ior of the condensate oscillations on and off resonance. 

In Section II we introduce briefly the kicked and 
double-kicked BEC systems. In Section III we intro- 
duce the time-dependent Bogoliubov method proposed 
by Castin and Dum and present numerics for the growth 
of non-condensate atoms. In Section IV we introduce a 
simple perturbative model, based on the one period time 
evolution operator for a kicked BEC. In Section V we 
show that the simple model and the time-dependent Bo- 
goliubov numerics give excellent agreement in the limit 
of weak kicks. In Section VI we consider the case T = 2tt 
with both numerics and the perturbative model and show 
that the instability border found in 0, @] is due to a novel 
type of compound Bogoliubov resonance. 



KICKED BEC SYSTEMS 

As in Q , we consider a BEC confined in a ring-shaped 
trap of radius R. We assume that the lateral dimension 
r of the trap is much smaller than its circumference, and 
thus we are dealing with an effectively ID system 



21] 



The dynamics of the condensate wavefunction at tem- 
peratures well below the transition temperature are then 
governed by the ID Gross-Pitaevskii (GP) Hamiltonian 
with an additional kicking potential: 



H = H GP + K cos f(t), 



where 



Hqp — — 



<:) 2 
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(1) 



(2) 



The short-range interactions between the atoms in 
the condensate are described by a mean-field term with 
strength g = 8N to tQ<sR/ r2 , where as is the s-wave scat- 
tering length, and N to t is the total number of atoms. For 



the QKR-BEC system, f(t) = J2„ S(t - nT ), whil e for 
the QKR2-BEC, 

/(*) = E P(* - nT ) - *(* - nT+ e )i > ( 3 ) 



where T is the total period of the driving; e <C T and 
thus the second kick nearly cancels the first. 

Experimental and theoretical studies of the double- 
kicked rotor have shown that its quantum behavior is 
largely determined by an effective kick strength K e = Ke, 
provided T ^> e. Here we take e = 1/25. Hence, while for 
the QKR-BEC, the value K = 1 represents a relatively 
large impulse for a kicked BEC, for a double kicked BEC, 
K = 1 in the numerics below corresponds to K e = 0.04, 
and represents only a very weak impulse. The reason for 
this is the near cancellation of consecutive kicks in each 
pair. 

This mechanism has certain analogies with the so- 
called "quantum antiresonance" investigated in for 
QKRs kicked at T = 2ir, consecutive kicks effectively 
cancel. This means that even large values of K ~ 1 and 
g > 1 represent only weak driving; for example, the in- 
stability border was found by Q to occur at g ~ 2 and 
K = 0.8. 



TIME-DEPENDENT BOGOLIUBOV METHOD 

The number of non-condensed atoms were calculated 
by making the usual Bogoliubov approximation, and fol- 
lowing the formalism of Castin and Dum [22| ■ This adap- 
tation of the Bogoliubov linearization for time-dependent 
potentials has been used in all studies to date of the dy- 
namical stability of kicked condensates 

HIE mill- The 

mean number of non-condensed atoms at zero tempera- 
ture is given by N ex (t) = Y!,kLi{ v k(t)\vk(t)) , where the 
amplitudes (v,k, Vfc) of the Bogoliubov quasiparticle oper- 
ators are governed by the coupled equations 



■t, d 
in— 

dt \v k 



H + gQ\i>\ 2 C 
-gQ*r 2 Q 



gQ^ 2 Q* 
H-gQ*\i>\ 2 L 



(4) 

In this expression, Q = I — \ip)(ip\ are projection op- 
erators that orthogonalize the quasiparticle modes with 
respect to the condensate 22j. We assume that at time 
t = 0, we have a homogeneous condensate i/jq = \/y/2~n. 
Further discussion of the theory is given in [23| . 

The regime of validity of the method is discussed in 
(22] ]. The method is valid in the weakly interacting limit 
1 » a? s p where p is the density. A limit is identified where 
this condition is satisfied, if one works with a constant 



totO'S 



thus the limit 



corresponds to N t , 



00. A further requirement is that condensate depletion 
remains negligible. This condition fails after a few kicks 
in exponentially unstable regions. Here the method is 
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employed only to identify the the parameter range for 
the onset of instability. We cut-off our calculations for 
N ex > 10 3 (a reasonable threshold for small depletion in 
a condensate with N tot ~ 10 5 ). 

In Figs[T] (a) and (b) we show the number of non- 
condensed atoms, N ex (t = NT), calculated from the 
Bogoliubov equations (0| after N = 200. For small 
K = 0.2, g = 1, a single resonance is seen at T ~ 10. 
For small K, resonances occur whenever the resonance 
condition 16] uj + toi = u\ ~ is satisfied, where 
n = 1, 2, 3.. is an integer and u>i is the eigenfrequency of 
the I — th collective mode. For larger K = 1, the figure 
shows that resonances are extremely dense and overlap 
with each other (and we show the behavior in this regime 
for T < 10). For overlapping resonances, unambiguous 
identification of each resonance is no longer possible. 
The key point here, however, is that in the stable regions 
outside the resonances, N ex remains very small even 
after prolonged kicking. 

FiglUb) shows oscillations of N ex , as a function 
of time, for weak K = 0.2, g = 1, close to the 
isolated resonance at T w 10. The three possi- 
ble regimes of: (non-resonant) weak quasi-periodic 
oscillations in time; (near-resonant) slower, large 
periodic oscillations; and (resonant) exponential 
growth are illustrated. The condensate energy, 

E(N) = fZ*Mr(N)(-$& + l\m)\ 2 )m) after N 
kicks, obtained from the GPE itself, is also shown, for 
comparison, in the inset: at resonance, large oscillations 
are also seen. 

Fig[2]shows the corresponding behavior for the double- 
kicked BEC, but now as a function of g, keeping T = 2, 
e = 1/25 constant and K — 1 or K = 5 (hence K e = 0.04 
or 0.2). The curve K e = 0.04 corresponds to weak im- 
pulses and shows two isolated Bogoliubov resonances. 
While values of g ~ 10 are large compared with current 
experimental values of g ~ 0.5 (see discussion of experi- 
mental g in [13| ) , resonances at small g ~ 1 more suitable 
for experimental spectroscopy can be excited by consid- 
ering larger T. The curve K e = 0.2 is in the overlapping 
resonance regime, so produces generalised instability. 

In order to understand the behavior at the resonances, 
we introduce in Section II a model for the time evolu- 
tion of perturbations from the kicked condensate, based 
on the usual linearization with respect to small pertur- 
bations. 
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FIG. 1: (a) Shows that for weak kicks (solid line), instability 
occurs only at one isolated Bogoliubov (1,1) resonance , where 
(n, I) denotes the n — th resonance of eigenmode I. "Up" ar- 
rows indicate onset of exponential instability; " down" arrows 
means stability is regained, g = 1. The total number of non- 
condensate atoms generated after 200 kicks, N ex (N = 200), 
is plotted as a function of kicking period T. For stronger 
kicks (dotted line; K = 1,T < 10) resonances proliferate and 
there is instability over almost all the parameter range, (b) 
Time-dependence near the(l, 1) resonance at T « 10 corre- 
sponding to Fig (a). Non-resonant (T = 13) curve shows weak 
quasi-periodic oscillations in N ex ; the near-resonant regime, 
T — 10 is characterized by slow, large oscillations; at reso- 
nance T = 10.5, there is exponential growth in N ex (t). Inset 
shows that the condensate energy (calculated from the GPE 
itself) has similar oscillations. 
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FIG. 2: double-kicked BEC (QKR2): Shows zones of insta- 
bility occur at Bogoliubov resonances. Condensate losses as 
a function of nonlinearity parameter g. "Up" arrows indicate 
onset of exponential instability; " down" arrows means sta- 
bility is regained. N ex (t = 1000) is plotted as a function of 
g (for T = 2, e = 1/25) for weak kicks (K = 1 so effective 
kick is K € = 0.04) and stronger kicks (K — 5 so effective kick 
R\ = 0.2). 



II: KICKED CONDENSATE MODEL 

The time-evolution of small perturbations of the con- 
densate itself arc described by an equation similar to 
EqHJ see 22]. We write the condensate wavefunction 
in the form ip = i/jq + 5ip, where i/jq is the unperturbed 



condensate and Sip represent the excited components. In- 
serting this form in the GPE and linearizing with respect 
to dtp, we can write: 
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where, 



C(t) 



H(t)+g\ip\ 2 



gr 2 

-H(t) - c 



(6) 



The analysis of condensate stability for a time-periodic 
system [l6[ reduces to the analysis of the operator C(t) 
over one period T. In general, for systems like BECs in 
modulated optical lattices, inter-mode coupling requires 
a detailed analysis of the instantaneous evolution. The 
nature of the (5-kicked potential permits considerable sim- 
plification. 

The effect of C(t) reduces to the free- ringing of the 
eigenmodes of the unperturbed condensate for period T, 
interspersed by instantaneous impulses which mix the 
modes. Even for an experiment (where the kicks are 
approximated pulses of very short, but finite duration) 
numerical time-propagation is avoided: intermode cou- 
pling occurs over a very short time-scale, during which 
eigenmode phases remain essentially constant. 

Excluding the kick term for the moment, we recall 
that the time propagation under Hqp can be analyzed 
in terms of the eigenmodes (uk(t), v k (t)) and eigenvalues 
of ui k (t) of the 2x2 matrix on the right hand side of 
Eq|6] Setting ip = l/s/2n, the matrix can be diagonal- 
ized and there are well-known analytical expressions for 
the unperturbed eigenmodes [l| 



(u k (t = 0),v k {t = 0)) 



ikl) 



/2tt 



(7) 



where U k + V k = A k , U k - V k = A k 1 , and A k = 

(^W + f)) 1/4 - 

In order to understand the behavior at the resonances, 
we introduce below a simple model using the eigenmodes 
EqJT] as a basis. Writing the small perturbation in this 
basis: 



ikO 



2tt 



-ikd 



(8) 

Neglecting the kick, evolving the modes from some ini- 
tial time to, each eigenmode (u k ,v k ) simply acquires a 
phase ie: 



b k (t)=b k (t )e-^ t - t °\ 



(9) 



where w k 



After a time interval T, a kick is applied which cou- 
ples the eigenmodes. Its effect is obtained by expressing 
the perturbation in a momentum basis, ip — ai(t)\l) 

where \l) — ^-7=, and we can restrict ourselves to the 
symmetric subspace a; — a_/ of the initial condensate 
(parity is conserved in our system) . Then, we can see by 
inspection that 



a k (t) = U k b k (t) + V k b*_ k (t). 



(10) 



Note that b k = b- k for this system. Conversely, the cor- 
responding amplitude b k in each eigenmode k is given 
trivially from EqJS] using orthonormality of the momen- 
tum states and the relation f7| — V k = 1, yielding 



b k (t) = U k a k (t)-V k a* k (t). 



(11) 



If the evolving condensate is given in the momentum 
basis, the effect of a kick operator U k i ck — e ±L a~ cose is 
well-known. The matrix elements: 



U n i = (n\U k ick\l) = Jn-i{K/Ti)i 

The J n -i are Bessel functions. 
The amplitudes ai (t) are given by 



±(l-7 



a n {t+) 



ai{t ), 



(12) 



(13) 



where a n (t + )/ai(t~) denotes the amplitude in state |n) 
just after/before the kick. 

We can now define a "time-evolution" operator 
£{T) = B^ 1 Cf ree (T)B Ukick-, where Cf ree denotes free 
ringing of the eigenmodes, B is the transformation from 
momentum basis to Bogoliubov basis and U k i ck is the ac- 
tion of the kick. A usual procedure for stability analysis 
of a driven condensate is to examine the eigenvalues of 
C'(T) to ascertain whether there is one (or more eigen- 
values) which have a real, positive component [16j |. ie 
whether they produce exponential growth in the ampli- 
tudes a±i. 

However, to compare with GPE numerics, we simply 
evolve the mode amplitudes in time over a few kicks and 
examine the overall condensate response to the kicking 
(in the limit of very weak kicking) . Hence we can evolve 
the amplitudes a; (t = NT) of the condensate perturba- 
tion from period N to period N + 1: 



&((N+ 1)T) = £{T) a(A0, 



(14) 



using only the simple analytical coefficients in Eg J13I 
and EqJHl provided we use the simple transformations in 
EasfTUland Eq[TT]to switch between the Bogoliubov mode 
basis and the momentum basis. C'(T) is non-unitary, but 
the method is quantitative in the perturbative limit pro- 
vided i\) ~ ipo, ie we assume ao(N) = ao(0) = 1. 

We calculate the average energy over the first few N 
kicks, (E(N)) = YltLi ^(*)- Slow, large amplitude os- 
cillations in E{t) yield a large (E(t)) and indicate a reso- 
nance. FigJ3^a) shows the QKR-BEC behavior, for equiv- 
alent parameters to Figflja). For low K — 0.2, there is 
the same single (1, 1) resonance at T sa 10 as in FigOTa). 
For higher K = 1 the method is far from quantitative: 
the model Eci U4l is only a valid means of time-evolving 
the perturbation over a few kicks for small K « 1 since 
it assumes the perturbed component is negligible; nev- 
ertheless, for K = 1 it illustrates the regime of dense, 
overlapping resonances. 
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In Fig|4|b) we compare the perturbative Eq. (fT4|) 
results with full GPE numerics for the first 20 kick 
pairs of the QKR2 in the limit of weak kicks. It shows 
remarkably good agreement. Moreover the scaling of 
the resonances with K is well described. The QKR2 
resonant Bogoliubov spectrum differs appreciably from 
the QKR case. FigfJJh) shows that for QKR2, even for 
low K = l,e = 1/25 ie K e w 0.04 and low g = 1, both 
I = 1 and I = 2 resonances are strongly excited. The 
QKR2-BEC resonance intensity depends strongly on K: 
the I = 2 resonances scale as K 4 , while the I = 1 scaling 
is closer to K 2 . In the full GPE numerics, the position 
of the maxima depends slightly on K and g 7 but remains 
within a few percent of the unperturbed value, even for 
longer kicking times if K e remains small. 

In the limit of weak driving, one can obtain explicit ex- 
pressions for the condensate wavefunction as a function of 
time. We assume that ao ~ \j\f2-K 3> a;^o- Then EqfT3l 
can be approximated by a n {t + ) w ai{t~) + Uio/(V2n). 
From Eg fTTI and Eq[5]we see that the amplitude accumu- 
lated over a single period in each eigenmode is 



h{N + 1) = h{N) + (UiUw - Wfo)e 



(15) 



Summing all contributions iteratively from t = 0, taking 
bi(0) = 0, we obtain 



b t {N) = (UiU l0 -ViUt ) 



n=N-l 
n=0 



(16) 



and so for ujiT w 2ir all the contributions add in phase, 
analogously to the well-known (but unrelated) resonances 
of the non-interacting limit [11| . 

We can write Yln=o e muJlT = e -^-l)<^$(*|^) 
where the $ function is: 



Nu) t T\ _ sm(NcJiT/2) 
~ sin (ujiT/2) 



(17) 



We thus expect oscillations in each set of ±1 momentum 
components of amplitude 



|2aKAT)| 2 oc4|^o| 2 $ 2 (^) 



(18) 



Off-resonance there will be quasi-periodic oscillations (in 
e.g. the condensate energy) from the superposition of 
contributions characterized by different eigenfrequencies 
u>l. Close to resonance, a single component dominates; if 
the I — th mode is resonant we can write w{T m 2nM + 25 
where 25 <C 1 is the de-phasing from resonance. Then 



\Uio[ 
5 2 



■sin 2 (iV(5), 



(19) 



and there are slow, periodic oscillations of large ampli- 

U w \" 



tude ~ 4 1 j° 1 , at a frequency 5 which is not related to 



any eigenmode frequency, but given rather by the de- 
phasing from resonance. 

The QKR2 resonant excitation spectrum is rather dif- 
ferent from the QKR, and is analysed further in the next 
section. 



Ill: RESONANCES OF THE QKR2-BEC 

In the limit K t — * 0, we can obtain analytical expres- 
sions for the BEC wavefunction of the double-kicked sys- 
tem. Firstly note that when ge -C 1, the non-linearity 
has little effect during the short time-interval e. Using 
the relation, 



= +i£cos0 -ip 2 th/2 -if-cos(e) 



r [p+K sin ( 



(20) 



the time evolution can be given as a 'one-kick' operator 



U(T)nUgl(T,0)e 



- [p-\-K sin f 



(21) 



In the limit pe ps 0, one can split the operators in Eq 1211 
and neglect a term K sin 9 ft to obtain the approximation 



U(T) 



. sin 2 6-iKeh cos $\ 



(22) 



leaving an effective single-kick quantum rotor with a kick- 
ing potential 



K 2 e 



■ sin 9 — iKeh cos 9 



J2S(t-NT). (23) 



JV 



The second term, curiously, appears as kicking potential 
with an imaginary, and h dependent, kick strength iKh. 
It is of purely quantum origin as it arises from the non- 
commutativity of p and sin#, i.e. 



iKh cos 9 — [K sm9,p\. 



(24) 



Nevertheless, as seen below, it is important for weak driv- 
ing as it controls the amplitude of the first excited mode 
l = ±l. 

The matrix elements of the modified kick Vkick, like 
those in Eq[13l are Bessel functions. Specifically, the 
effect of Vkick on the condensate amplitudes a/ is given 
by 



*n(* + ) = X! Unl a '(* )' 



(25) 



where U n i = X) m i" 1 m Jm(^ff) J n -l-2m(iKe), and 
&n(t ) indicates momentum amplitudes before(-) and 
after(+) the kick, as in Eq[l3] Since Ke <C 1 and 
J\ n \>i( z ) — 0, only Bessel functions of low order (m = 
or 1) will be non-negligible, and we can use the small- 
argument approximations for them, namely Jo(z) ~ 1, 
J±i(z) w ±z/2. 
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FIG. 3: Average energy (E) after 40 kicks. The dashed lines 
indicate the model of Eq |14l all other plots use full numerics. 
The label (n, I) denotes n — th resonance of mode I. Reso- 
nances of the QKR-BEC for parameters comparable to Fig. la. 
For low g = 1, K = 0.2, only the single isolated (1,1) reso- 
nance is seen. For higher K = 1, resonances proliferate and 
overlap. 
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FIG. 4: Comparison between full GPE numerics and the 
model of Eg 1141 for the QKR2-BEC, showing excellent agree- 
ment. Average energy (E) after 20 kick-pairs. The label (n, I) 
denotes n — th resonance of mode I. g = 1 and e — 1/25 so 
K = 1 corresponds to effective kick strength K t = 0.04. For 
low K, the I = 1 resonance amplitudes scale as ~ K 2 while 
those of the I = 2 modes scale as ~ K 4 . 



resonance of either the I = 1 or I = 2 modes, the regime 
illustrated in Fi^b). 

Similarly as for the QKR-BEC, we can sum all the con- 
tributions to obtain an approximate analytical expression 
for the evolving condensate wavcfunction including ex- 
cited modes I = ±1 and I — ±2, 

*P(N) « -^1 + Ci^cos0 + C 2 ; ^cos20]. (28) 

Z7T Z "til 

where 

d = -$(A r o7i)[cos(^- l)o7i -^ 2 sin(Ar- l)oJi], 
C 2 = ${N£i)[J$wn.{N - l)cJ 2 + icos(N- l)cJ 2 , 

and <Jjj = ujjT/2. 

Eq r251 shows that the amplitudes | &i | 2 and |a 2 | 2 scale 
as K 2 and K 4 respectively, as seen in the numerics in 
FigfJJb). FigEI a ) shows that Eqf25]gives excellent agree- 
ment with GPE numerics, giving accurately the non- 
resonant quasi-periodic condensate oscillations. Near the 
/ = 2 resonance of Fig(2l FigO^b) confirms the QKR2 
condensate oscillations (obtained from the GPE) scale 
quite accurately as oc sin 2 NS as expected from Eq[T2] 
and Eq[28l 

FigEJc) shows that, near-resonance, there are corre- 
sponding large oscillations in the non-condensate num- 
bers calculated from EqHl Near-resonance, N ex increases 
quadratically with time, on-resonance, the increase is ex- 
ponential. 

IV: BOGOLIUBOV RESONANCES FOR T = 2tt 



Then, if the condensate is relatively unperturbed, the 
main effect of the kick will be to simply excite a small 
amount of I = ±1 and I = ±2 from the |0) state 



e -%v k , ck ^ w e -±;V kick = ^ Ui \0) 

I 

where 



(26) 



1 



, iKe\ , v /K 2 e . 



2tt 4 1 1 8h 1 1 



(27) 



We obtain a similar equation to the QKR-BEC for the 

mode amplitudes, i.e. 

6, (AO = (UiUio ~ Wf ) exp[inuiT\. 

But if only the lowest excited modes are significant, 
then, in particular, 



h(N) = -^(U 1 -V 1 ) EnJo 1 expIinwrT] and 
b 2 (N) = i ^(U 2 + V 2 ) En=o exp[mw 2 T]. For lo{F « 2tt 
all the contributions add in phase and we will have a 



The kick period T — 2it, in a non-interacting system 
of cold atoms (i.e. g = 0) corresponds to a so-called 
"quantum anti-resonance" where the cold atom cloud 
exhibits periodic (period-2) oscillations. Hence the 
isolated Bogoliubov resonance regime to higher K than 
would be expected for generic T. The effect of a non-zero 
g for T = 2ir was investigated in t 6|. An instability 
border occurring at a critical value of nonlinearity, e.g. 
for g ~ 2 at K — 0.8, was identified where the growth on 
non-condensate particles with time became exponential. 

In Figl^a) we investigate the behavior near critical 
g, for K = 0.8. We see that if a wider range of g is 
considered, the stability border is also a resonance: the 
condensate rapidly recovers stability after the instability 
border is passed. The condensate is exponentially unsta- 
ble for g ~ 2 — > 2.6, but is quite stable for both g = 1.5 
and g = 3, as shown. FigO^b) shows oscillations in the 
condensate energy, as a function of time; a smoothed plot 
is also shown. For g — 1.5 and g = 2.8 (off-resonance) 
the smoothed plots are flat; for g — 2.2 and g = 2.5 
(near-resonant), slow deep oscillations are apparent. 
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Off-resonance 




20 40 60 
near resonance 



80 N 100 









200 time (N) 400 



FIG. 5: Test of perturbative model, (a) Condensate energy 
oscillations from GPE numerics and EqfSU K € — 0.04, g — 2, 
T = 2. Beating between modes 1 and 2 is very accurately de- 
scribed by Eq |28l (b) Behavior of I — 2 resonance of Fig. 1(b) 
K t = 0.04,T = 2 and g ps 9.5. As the resonance is approached 
the amplitude of the oscillations is proportional to the square 
of their wavelength, i.e. E(t = NT) cx K i -^ 7 sin 2 N8 where 
28 is the distance from the resonance peak, (c) Correspond- 
ing number of non-condensate atoms from Eq|4] 



The behavior is analogous to that of generic T; how- 
ever, the analysis of the condensate resonances for T = 
2ir is less straightforward: the strongest resonances, even 
for low K < 2, do not in fact occur for uo{T as 2irM, 
where M = 1,2,3.... 

A significant difference between generic T and T = 2n 
is that, for the generic case, if we write 



10 
2 




i . i . i . i . i 

0.5 1 1.5 2 2.5 3 
g 



50 100 

Time(kicks) 



50 100 
Time (kicks) 



FIG. 6: (a) Non-condensate particles for kicking period T = 
2tt, K — 0.8, g ~ 2. The inset shows the rate of exponen- 
tial growth of non-condensate atoms; zero denotes polynomial 
growth or less. The graph shows this instability border is a 
resonance: the condensate is unstable for g — 2 — 2.5 but is 
stable for g = 1.5 and g = 3.0. (b) Energy oscillations as a 
function of time; smoothed plots are also shown. Before and 
after the resonance (g = 1.5 and g = 2.8) the smoothed plots 
are flat. Near-resonance, (g = 2.2 and g = 2.5) the energy 
shows the characteristic slow, deep resonant oscillations. 



ujiT as 2-kMx + 25(1) 



(29) 



we see that for arbitrary generic T, the distance from the 
nearest resonance, for the different modes, depends on I. 
In contrast, for T = 2ir, for large I (i.e. I > 3) we find 
ujiT w (I 2 + £)n; in other words, the de-phasing from 
the nearest resonance (and hence the period of the mode 
oscillations) is similar (either 2<5 w ^ or 25 w 1 — f ) 
for for all modes. So all mode oscillations for high / are 
approximately in phase with each other. 

For K = 0.8, only low modes I = 1, 2 are significantly 
populated. These low modes (/ = 1 and I — 2) are only 
in phase with each other at certain precise values of g, T. 
For these parameters, the model of EafT3l predicts large 
resonances whenever the condition (cji +ui2)T w 2-7rM is 
satisfied. In particular, for the resonance near g as 2, we 
find that for the I = 1 mode, oj\T as (1 — 25)2tt while for 
the I = 2 mode u 2 T w (2 + 25)2tt, with 25 w .25. 



These results suggest that "two-mode resonances" , i.e. 
synchronized oscillations of pairs of the lowest excited 
modes are the dominant mechanism for T — 2ir (NB 
this could be viewed as a "three-mode" resonance, if we 
include the lowest, initial mode, but ujq = for our sys- 
tem). They account for the shifting position of the crit- 
ical instability border found by in the T = 2tt case. 
For example, for slightly higher kick strengths, such as 
K ~ 2, a resonance appears for g as 1.65 corresponding 
to (k>2 + ^z)T as 2nM, which accounts for the displace- 
ment of the instability border to lower values of g. Note 
that the resonance positions in the full numerics are in- 
dependent, whereas in the perturbative model of EafLlI 
this dependence is neglected; the model is only valid for 
very small K. 
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V: CONCLUSION 

In conclusion, we have shown that exponential insta- 
bility in kicked BECs is related to parametric resonances, 
ie driving of low-lying collective modes at their natural 
frequencies, rather than to chaos in the underlying mean- 
field dynamics [24l |. 

The signature of this process is in the onset of slow, 
large amplitude periodic oscillations in the condensate 
energy as well as the number of non-condensate atoms 
calculated from the time-dependent Bogoliubov formal- 
ism, as a resonance is approached. The resonances pro- 
liferate and overlap for large kick-strengths K, leading 
to instability over wider ranges of K and g. The time- 
dependent Bogoliubov approximation used here and in 
all other previous studies is only valid in regimes where 
the condensate depletion is negligible; for realistic con- 
densates analysis of the dynamics in the narrow (for 
weak driving) windows of parametric instability, would 
require other approaches beyond Bogoliubov. However, 
away from these windows, the kicked condensate remains 
stable and relatively unperturbed, even after prolongued 
kicking. 
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